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LONG-TERM  GOALS 

The  long-term  goal  is  to  implement  input  and  dissipation  source  functions,  based  on  advanced 
understanding  of  physics  of  air-sea  interactions,  wave  breaking  and  swell  attenuation,  in  wave- 
forecast  models. 

OBJECTIVES 

The  objectives  are  to  use  new  observation-based  source  terms  for  the  wind  input,  wave-breaking 
(whitecapping)  dissipation  and  swell  decay  in  the  third-generation  models  WAVEWATCHIII  and 
SWAN.  Calibration  and  performance  of  the  source  functions  have  to  satisfy  a  set  of  physical 
constraints,  and  methodology  is  to  be  developed  to  enable  testing  the  source  functions  separately 
before  they  are  blended  in  the  full  model.  Verification  is  to  be  conducted  by  means  of  academic 
tests  and  hindcasting  real-life  scenarios  defined  by  the  end  users  from  the  US  Navy,  Army  and 
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NOAA,  to  include  deep  and  finite-depth  conditions,  closed  seas  (no  swell)  and  open  ocean,  extreme 
weather  events,  and  global  simulations. 

APPROACH 

Physics  of  two  primary  source/sink  terms  employed  by  the  operational  models,  namely  wave¬ 
breaking  energy  dissipation  and  wind-to-wave  energy  input  have  not  been  updated  for  decades.  In 
the  meantime,  the  new  physics  is  available.  For  the  first  time  under  field  conditions,  in  the  course 
of  the  ONR  Lake  George  (Australia)  project,  estimates  of  the  spectral  distribution  of  the  wave¬ 
breaking  dissipation  were  obtained,  and  measurements  of  the  wind-input  spectral  function  were 
conducted  at  moderate-to-strong  wind  forcing  (Young  et  al.  2005).  Corresponding  outcomes  were 
parameterised  as  source  functions  suitable  for  spectral  wave  models,  and  both  exhibit  a  number  of 
physical  features  presently  not  accounted  for. 

For  the  dissipation,  these  are  threshold  behaviour  of  breaking/dissipation  in  terms  of  wave- 
steepness/spectral-density,  cumulative  effect  at  scales  smaller  than  the  spectral  peak,  and  direct 
coupling  of  the  dissipation  rates  with  input  rates  at  very-strong/extreme  wind  forcing.  Bi-modal 
directional  distribution  of  the  dissipation  was  also  observed  (Babanin  et  al.  2001,  Babanin  and 
Young  2005,  Young  and  Babanin  2006a,  Babanin  et  al.  2010).  None  of  these  features  are  present  in 
the  currently  employed  dissipation  functions  used  for  operational  wave  forecasting. 

The  new  wind-input  features  are  non-linear  behaviour  of  the  input  term  (that  is,  the  input  rates 
depend  on  wave  steepness)  and  full  flow  separation  in  extreme  conditions  (that  is,  relative  slowing 
down  of  the  wind-wave  exchange  in  steep- waves/ strong- winds  circumstances)  (Donelan  et  al. 

2005,  2006,  Tsagareli  et  al.  2010).  Enhancement  of  the  wind  input  due  to  wave  breaking  was  also 
observed,  quantified  and  parameterised  (Babanin  et  al.  2007a). 

While  the  wind  input  and  dissipation  are  the  main  source/sink  energy  terms  in  the  model,  the  latter 
has  to  be  subdivided  into  separate  terms:  one  term  in  case  of  breaking  waves  (for  wind-generated 
waves),  and  another  for  non-breaking  waves  (swell).  The  first  one  (whitecapping  dissipation)  turns 
zero  once  the  spectral  density  at  a  particular  scale  drops  below  the  threshold  as  mentioned  above, 
and  a  much  weaker  dissipation  accompanies  swell  propagation  across  the  ocean.  Recently,  two  new 
swell-dissipation  terms  were  suggested,  which  imply  different  physics:  interaction  of  swell  with 
atmospheric  boundary  layer  (Ardhuin  et  al.  2009)  and  with  background  turbulence  (Babanin  2006, 
2011a). 

The  project  will  use  both  the  new  source  terms  and  a  new  approach  to  their  validation.  The  main 
feature  of  the  approach  is  stringent  physical  constraints  on  the  momentum/energy  fluxes  in  and  out 
the  wave  system:  that  is,  the  integrated  momentum  wind  input  must  not  exceed  independently 
known  total  stress,  and  the  integrated  dissipation  must  constitute  the  experimentally  known 
proportion  with  respect  to  the  total  input.  For  the  total  dissipation,  independent  parameterisations 
based  on  profile  measurements  of  volumetric  dissipation  rates  are  also  available.  Such  constraints, 
first  of  all,  are  necessary  to  make  the  source  functions  physically  consistent,  and  most  importantly, 
they  allow  to  calibrate  the  input  and  dissipation  one  by  one,  before  they  incorporated  in  the  model 
where  their  contributions  are  not  possible  to  separate  (Babanin  et  al.  2005,  2010,  Tsagareli  et  al 
2010,  Rogers  et  al.  2011). 

Other  constraints  were  also  formulated  and  calibration  dependences  for  initial  academic  tests  were 
selected  (see  Work  Completed  below).  For  practical  testing  and  hindcasting,  a  set  of  field  sites  and 
datasets  were  chosen  which  include  Lake  Michigan  (deep  water,  no  swell,  Rogers  et  al.  2011), 

Lake  George  (finite  depth,  no  swell,  Young  and  Verhagen  1996a, b,  Young  et  al.  1996,  Young  and 
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Babanin  2006b,  2009,  Babanin  and  Makin  2008,  Young  2010),  Gulf  of  Mexico  (open  sea,  deep-to- 
finite  depths,  swell,  hurricanes,  IPET  2006,  Smith  2006)  and  a  selection  of  tropical  cyclones  from 
the  Australian  region  (Young  2006,  Babanin  et  al.  201  la).  Measurement  data  are  available  for  all 
the  chosen  sites.  Once  the  model  is  fully  tested,  global  simulations  for  selected  years  will  be 
performed  and  compared  with  the  altimeter  database  (Zieger  et  al.  2009,  Vinoth  and  Young  2011, 
Young  et  al.  201 1)  and  with  NOAA  simulations  based  on  the  previous  version  of  WAYEWATCH 
III. 

The  research  group  includes  academics  from  Swinburne  University  of  Technology  (SUT),  US 
Naval  Research  Laboratory  (NRL),  US  Army  Corps  of  Engineers  (CHL)  and  US  Weather  Service 
(NOAA).  SUT  group  consists  of  Alex  Babanin  and  Ian  Young  who  had  developed  the  new 
dissipation  term  and  participated  in  development  and  testing  of  the  new  wind-input  term  in  the 
course  of  the  earlier  ONR  project,  and  Stefan  Zieger,  Research  Fellow  employed  on  the  current 
project.  This  group  will  conduct  bulk  of  the  academic  research  and  implementation  of  the  new 
source  terms  into  WW3,  the  latter  in  collaboration  with  Hendrik  Tolman  from  NOAA  and  Erick 
Rogers  from  NRL.  Erick  Rogers,  and  Jane  Smith  from  CHL  will  contribute  to  the  project  goals  in 
the  development  and  verification  of  the  new  physical  formulations,  to  replace  the  formulations 
currently  used  in  Navy  and  Army  operational  models.  A  specific  goal  of  NRL  and  CHL  is  to  create 
models  that  are  physically  consistent  with  what  we  know  about  the  real  ocean,  while  at  the  same 
time  ensuring  that  the  new  models  are  optimal  for  Navy  and  Army  applications.  NRL  and  CHL  will 
implement  the  new  input  and  dissipation  terms  into  SWAN,  apart  from  WW3,  and  will  participate 
in  the  development  of  validation/calibration  cases,  and  in  particular  those  of  extreme  conditions 
observed  in  the  Gulf  of  Mexico. 

WORK  COMPLETED 

This  is  a  report  for  the  second  year  of  the  project.  In  the  first  year,  the  formulations  for  the 
observation-based  source  terms,  were  validated,  individually  calibrated,  and  tested  together  in  a 
two-dimensional  wave  research  model  with  exact  computations  of  the  non-linear  interaction  term 
(Tsagareli  et  al.  2010,  Babanin  et  al.  2010).  The  new  source  terms  and  the  physical-constraint 
approach  were  implemented  in  SWAN  (Rogers  et  al.  2011)  and  tested  by  means  of  the  Lake 
Michigan  and  a  selection  of  the  Gulf  of  Mexico  cases.  Field  data  sets  for  further  model-testing  were 
selected  and  prepared:  Lake  George;  Black  Sea;  global  altimeter  data  base. 

During  the  second  year,  the  main  aim  of  the  project  was  implementation  of  the  source  functions 
into  WW3.  This  is  done,  and  academic  testing  and  initial  field  validation  of  the  model  were 
concluded.  In  addition,  the  model  performs  an  automatic  self-correction  routine  by  comparing  the 
input  total  stress  with  the  integral  of  the  wind  input  over  the  computed  spectrum  at  each  time  step. 

The  academic  tests  covered  ability  of  the  new  version  of  the  model  to  reproduce  observational 
deep-water  dependences  for  integral  parameters,  for  spectral  parameters  and  for  directional 
parameters  of  wave  fields  (Babanin  and  Soloviev  1998a, b).  New  tests,  intended  on  separate 
validation  of  the  wind  input  function  and  the  whitecapping  dissipation  term  were  developed.  The 
field  validation  of  the  model  was  conducted  by  means  of  hindcasting  for  a  selection  of  Lake 
Michigan  storms  (wind  and  wave  data  are  available,  no  swell  conditions). 

Research  on  the  other  topics  of  the  project  continued.  On  wave  breaking  and  dissipation,  book  was 
published  (Babanin  2011a),  breaking  in  directional  wave  fields  (Babanin  et  al.  2011b,  Toffoli  et  al. 
201  la),  breaking  in  spectral  environments  (Chalikov  and  Babanin  2011)  and  wind  influence  on  the 
severity  of  the  breaking  (Galchenko  et  al.  2011)  were  investigated.  The  wind- input  constraint,  sea 
drag,  and  its  dependence  on  the  wave  breaking  (Babanin  201  la),  on  the  directional  spreading  of 
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waves  (Ting  et  al.  2011)  and  on  wind  gustiness  and  air  humidity  (Le  Roy  2011)  were  studied.  New 
dynamic  wave -bottom  friction  routine  was  suggested  and  tested  in  SWAN,  which  accounts  for 
ripple  formation  and  dissolution,  grain  size  of  the  sediment  (Smith,  G.  et  al.,  2011).  Continued  was 
research  of  the  wave  induced  turbulence  and  its  role  in  the  swell  dissipation  (Babanin  2011a, 
Pleskachevski  et  al.  2011). 

Particular  attention  was  paid  to  preparing  for  future  hurricane  modelling  and  global  simulations. 
Sensitivity  tests  for  the  source  functions  in  extreme  conditions,  including  bottom-limited  breaking, 
diffraction,  refraction  and  numerical  schemes  were  conducted  by  means  of  simulating  Typhoon 
Krosa  (Babanin  et  al.  2011a).  New  technique  for  measuring  spray  concentrations  in  the  course  of 
the  planned  Australian  tropical  cyclone  field  experiment  was  developed  (Toffoli  et  al.  201  lb). 
Change  of  dynamical  regime  of  air-sea  interactions  at  hurricane-wind  threshold  is  described 
(Babanin  201  lb).  The  global  altimeter  database  of  Zieger  et  al.  (2009)  was  used  for  investigating 
the  global  trends  of  waves  and  winds  (Young  et  al.  2011)  and  for  estimating  their  extremes  (Vinoth 
and  Young  2011),  over  the  past  25  years. 

RESULTS 


The  results  already  submitted,  accepted  for  publication  or  published  and  their  significance  are 
outlined  in  sections  Work  Completed,  Impact/ Applications  and  Publications  of  this  report.  In  this 
section,  due  to  limited  space,  we  will  indicate  on  results  achieved,  but  not  published  yet,  and  will 
highlight  most  relevant  published  results. 


The  academic  evolution  tests  of  new  wind-input  and  breaking-dissipation  functions  were  done  by 
two  means,  with  DIA  parameterisation  of  the  nonlinear  interactions  and  with  exact  computations  of 
these  interactions  by  means  of  WRT  code  (Tracy  and  Resio  1982).  The  latter  provides  accurate 
estimates  of  energy  fluxes  within  the  wave  system,  in  addition  to  those  due  to  wind  and  breaking, 
but  is  not  possible  to  employ  operationally,  and  the  former  is  approximate,  but  is  routinely 
employed  in  the  operational  forecast.  Results  are  summarised  in  captions  to  Figs.  1  and  2,  and  Fig. 

3  demonstrates  application  of  these  functions  to  Lake  Michigan  and  comparisons  with  buoy 
measurements. 


Fig.  1.  Total  wave  energy  (dimensionless)  versus  dimensionless  peak  frequency,  for  different 
source  functions  as  indicated  in  the  legends.  The  dependences  are  normalised  by  the 
observational  parameterisation  of  Babanin  and  Soloviev  (1998a),  which  therefore  equals  unity  in 
Figure,  with  shaded  areas  showing  95%  confidence  limits.  Green  lines  indicate  Tolman  and 
Chalikov  (1996)  source  functions,  blue  lines  WAM3  and  red  lines  WAM4  physics,  respectively. 
Suit  of  black  lines  correspond  to  Donelan  et  al.  (2006)  wind  input  in  the  Tsagareli  et  al.  (2010) 
version  and  to  the  four  versions  of  Babanin  and  Young  (2005)  dissipation,  following  Rogers  et 
al.  (2011).  UL4M4  is  the  preferred  version  for  the  subsequent  modelling.  (Left)  Exact  nonlinear 
term  is  used.  (Right)  DIA  nonlinear-interaction  term  is  used 
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The  integral  growth  curves  are  usually  subject  to  main  scrutiny  and  tuning  and  all  the  source 
functions  perform  reasonably  well  with  respect  to  the  experimental  dependence.  Since  the  tuning  is 
typically  done  with  the  use  of  DIA  in  evolution  simulations,  it  is  noticeable  in  Fig.  1  that 
performance  of  some  source  terms  actually  deteriorates  once  the  exact  nonlinear  term  is  employed. 

In  Fig.  2,  evolution  results  are  shown  with  the  use  of  WRT  nonlinear  term  only.  For  the  level  of  the 
spectrum  tail  (left),  the  new  dissipation  functions  lead  to  overestimating  this  level  by  comparison 
with  measurements  and  to  corresponding  bias  of  high  moments  in  the  hindcast  of  Lake  Michigan 
storms  in  Fig.  3.  In  Fig.  2  right,  Donelan  et  al.  (2006)  input  (top)  and  UL4M4  dissipation  (bottom) 
are  compared  with  the  same  set  of  source  functions  as  before.  The  differences  are  essential,  from 
overestimating  to  underestimating  by  a  factor  of  2  or  more. 


VMM 3  XNL 

—  —  TC96  XNL 

—  —  VMM 4  XNL 

-  STXXNLDL1M1 

-  SIX  XNLUL1M4 

.  STX  XNI  IJI  7  M? 

—  -  STX  XNL  UL4M4 


dimensionless  frequency  f.i( 


Fig.  2.  Evolution  tests  with  WRT  nonlinear  term  and  the  same  selection  of  other  source  functions 
as  in  Fig.l.  (left)  Spectrum  tail  level  versus  dimensionless  frequency.  The  dependences  are 
normalised  by  observational  parameterisation  of  Babanin  and  Soloviev  (1998a),  which  therefore 
equals  unity  in  Figure,  with  shaded  areas  showing  95%  confidence  limits,  (right)  Comparison  of 
wind  input  (top)  and  dissipation  (bottom)  functions,  normalised  by  Donelan  et  al.  (2006)  input 
function  and  UL4M4  dissipation  function,  respectively.  Uw  /  cp  — 1.9,  bottom  scale  is  relative 

frequency  f/fp. 


In  Fig.  3,  the  new  source  functions  are  used  for  a  Lake  Michigan  hindcast.  Since  the  integral 
dependences  are  reproduced  well  in  academic  tests  (Fig.  1),  so  are  the  integral  wave  properties  in 
Fig.  3  (wave  height,  peak  period,  mean  direction).  Since  the  tail  is  overestimated  in  academic  tests 
in  Fig.  2,  the  weighted  periods  are  expectedly  biased  in  Fig.  3.  It  is  now  known  experimentally  that 
at  strong  forcing  at  a  particular  scale  the  extra  wind  energy  flux  cannot  be  effectively  transferred  by 
nonlinear  interactions  and  a  part  of  it  is  dissipated  locally  (Babanin  and  Young  2005,  Babanin  et  al 
2007b).  This  effect  should  reduce  the  tail  level  and  is  now  parameterised  in  order  to  optimise  the 
tail  and  high-moment  behaviour  (not  shown). 
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Lake  Michigan  -  STX4  DIA 
NDBC  directional  wave  spectra 


NDBC  directional  wave  spectra 


Fig.3.  Comparison  between  the  model  (Donelan  et  al.  (2006)  input  and  UL4M4  dissipation  of 
Rogers  et  al.  (2011))  and  NODC  buoy  45007 from  a  3  day  hindcast  (12-15  October  2002).  DIA 
nonlinear  term  is  used.  Legends  in  the  panels  show  parameters  compared  and  statistics  of  the 
comparisons.  (Top  row,  left  to  right)  Wave  height,  Mean  wave  period,  Mean  wave  direction, 
respectively.  (Bottom  row,  left  to  right)  Peak  period,  T02  period,  7I0  period,  respectively 

One  of  major  developments  over  the  year  was  use  of  the  Swinburne  University  altimeter  database 
(Zieger  et  al.  2009)  for  investigating  changes  of  mean  and  extreme  values  of  winds  and  waves  since 
1985,  which  will  then  be  employed  in  global  testing  of  the  new  version  of  the  model.  Fig.  4, 
reproduced  from  Science  paper  by  Young  et  al.  (201 1)  shows  growth  of  their  values  globally. 
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99th  percentile  wind  speed  {1991-2008} 


Fig.4.  Colour  contour  plots  of  the  99th  -percentile  trend  (per  cent  per  year,  1985-2008  period, 
colour  scale  is  shown  on  the  right).  Points  that  are  statistically  significant  according  to  Seasonal 
Kendall  test  are  shown  with  dots,  (top)  Wind  speed,  (bottom)  Wave  height 

IMPACT/APPLICATIONS 

Paper  by  Rogers  et  al.  (2011)  is  implementation  of  the  new  wind-input  and  whitecapping- 
dissipation  functions  in  SWAN  model.  These  are  functions  based  on  observations  and  advanced 
understanding  of  respective  physics,  which  have  a  potential  to  become  operational  for  wave- 
forecast  models.  The  paper  also  suggests  and  employs  new  approach  to  verifying  the  source 
functions  individually  at  each  computation  step,  which  approach  has  a  potential  to  replace  the 
existing  tradition  of  testing  the  models  based  on  bulk  performance  of  the  energy  sources  together. 

New  swell-decay  function  developed  by  Babanin  (201  la)  is  based  on  measurements  and  modelling 
of  wave -turbulence  interactions,  a  phenomenon  not  accounted  for  before,  ft  provides  estimates  of 
swell  dissipation  consistent  with  field  observations  of  Ardhuin  et  al.  (2009)  and  is  a  critical  new 
term  for  the  models  with  threshold  behaviour  of  whitecapping  dissipation  where  such  dissipation 
ceases  once  the  wave  spectrum  is  below  the  threshold.  Dissipation  due  the  wave-induced 
turbulence  is  most  important  in  the  context  the  upper-ocean  mixing  and  sediment  suspension,  and 
this  topic  is  further  investigated  by  Pleskachevsky  et  al.  (2011). 

New  bottom-friction  routine  was  developed  and  introduced  in  SWAN  which  accounts  for  sediment 
size  and  provides  dynamic  coupling  of  the  ripple-caused  friction  with  the  waves  (Smith,  G.  et  al. 
2011).  Simulations  of  Typhoon  Krosa  by  means  of  SWAN  and  WWM11  revealed  bottom-limited- 
breaking  terms  as  most  problematic  in  finite  depths  (Babanin  et  al.  2011a).  These  results  have 
essential  implications  for  finite-depth  modelling. 
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Topic  of  extreme  wind- wave  interactions  is  most  important  for  the  project  in  general,  as  major 
hindcast  testing  cases  will  be  those  for  hurricanes  in  North  America  and  tropical  cyclones  in 
Australia.  Babanin  (201  lb)  analysed  this  topic  in  general  and  found  transition  to  a  new  dynamics  in 
all  three  relevant  media:  in  the  atmospheric  boundary  layer,  on  the  ocean  surface,  and  below  the 
surface,  at  the  hurricane-threshold  winds  of  30-35m/s.  Toffoli  et  al.  (2011b)  developed  a  new 
method  of  measuring  spray  in  field  experiments  in  hurricanes,  which  spray  is  a  critical  feature  for 
hurricane  dynamics  according  to  many  theoretical  models. 

Book  of  Babanin  (201  la)  published  by  Cambridge  University  Press  is  the  first  book  on  the  wave 
breaking  and  dissipation.  It  also  analyses  features  and  phenomena  in  the  atmospheric  boundary 
layer  and  below  the  ocean  surface  related  to  the  breaking,  including  those  phenomena  and 
transitions  which  occur  in  extreme  weather  conditions. 

Publications  of  Babanin  et  al.  (201  lb, c),  Chalikov  and  Babanin  (2011),  Galchenko  et  al.  (2011), 
Toffoli  et  al  (201  lb)  further  develop  the  topic  of  wave  breaking.  Quantitative  criteria  for  wave 
breaking  in  directional  fields  and  for  wave  trains  with  continuous  spectrum  are  analysed,  influences 
of  the  wind  forcing  on  wave  breaking  severity  are  studied.  Potential  impact  of  this  research  is 
significant  across  variety  of  related  topics,  including  modelling  the  wave  dissipation. 

The  sea-drag  was  shown  to  be  affected  by  wave  breaking  (Babanin,  2011a)  and  to  depend  strongly 
on  the  directional  spreading  (Ting  et  al.  2011),  on  the  gustiness  and  humidity  (Le  Roy  2011).  In  the 
context  of  the  present  project,  the  drag  is  used  as  a  wind-input  constraint,  but  impact  of  potentially 
improved  sea-drag  parameterisations  is  much  broader  in  the  air-sea  interaction  modelling  of  all 
scales. 

Essential  outcomes  came  from  analysis  of  altimeter  satellite  database  of  Zieger  et  al.  (2009).  Young 
et  al.  (2011)  revealed  global  positive  trends  for  both  mean  and  extreme  values  of  surface  winds  and 
wave  heights,  and  Vinoth  and  Young  (2011)  provided  estimates  of  global  extremes  for  such  winds 
and  waves.  These  results  will  be  used  for  global  testing  of  the  model  with  new  source  functions  at 
the  final  stage  of  the  project. 

Finally,  but  not  the  least  important,  is  impact  of  the  project  on  research  training.  4  PhD  and  1 
Master  theses  have  been  defended  or  submitted  on  the  topics  of  wave  breaking,  sea  drag,  altimeter 
data,  wave-bottom  interaction,  whose  results  were  influenced  and  were  used  in  the  present  project. 

TRANSITIONS 

The  input  function  (Tsagareli  et  al.  2010)  and  dissipation  (Babanin  et  al.  2010)  are  used  in  research 
wave  models  at  the  University  of  Darmstatd,  Germany  and  the  National  Cheng  Kung  University, 
Taiwan 

RELATED  PROJECTS 

Ardhuin  et  al.  “Ocean  Wave  Dissipation  and  Energy  Balance  toward  Reliable  Spectra  and  First 
Breaking  Statistics”.  NOPP  project,  implements  new  dissipation  function  based  on  similar 
physical  principles.  Joint  publications  (Ardhuin  et  al.,  2010,  Filipot  et  al.,  2010) 

Babanin,  A.V.,  Phillips,  W.R.C.,  Ganopolski,  A.  “Wave-induced  upper-ocean  mixing”,  Australian 
Research  Council  (ARC)  Discovery  grant.  Investigation  of  the  wave-induced  turbulence  and 
non-breaking  dissipation  (Babanin  2011,  Pleskachevsky  et  al.  2011) 
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Babanin,  A.V.,  Young,  I.R.,  Phillips,  W.R.C.,  Donelan,  A.M.,  Makin,  V.,  Ardhuin,  F.  “Oceanic 
Conditions  within  Extreme  Tropical  Cyclones”,  ARC  Linkage  grant.  Field  investigation  of 
wind  input  and  dissipation  in  extreme  conditions  (Toffoli  et  al.  201  lb) 

Young,  I.R.,  Babanin,  A.V.  “A  Global  Satellite  Altimeter  Database  for  Ocean  Engineering 

Applications”,  ARC  Linkage.  Creation  of  the  global  altimeter  wind  and  wave  database  over 
the  period  1985-2008  and  investigation  global  and  regional  trends  (Vinoth  and  Young  2011, 
Young  et  al.  2011) 

Young,  I.R.,  Babanin,  A.V.,  Stiassnie,  M.A.,  Greenslade,  D.J.  “Numerical  modelling  of  extreme 
waves  generated  by  tropical  cyclones”,  ARC  Discovery.  Modelling  tropical  cyclones, 
investigation  of  the  nonlinear  source  term  for  spectral  models  (Babanin  et  al.,  2011a) 
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